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Abstract 

We present holographic computations of the time-dependent chiral magnetic 
conductivity in the framework of gauge/gravity correspondence. Chiral magnetic 
effect is a phenomenon where an electromagnetic current parallel to an applied mag- 
netic field is induced in the presence of a finite axial chemical potential. Motivated 
by a recent weak-coupling perturbative QCD calculation, our aim is to provide a 
couple of complementary computations for strongly coupled regime which might be 
relevant for strongly coupled RHIC plasma. We take two prototypical holographic 
set-ups for computing chiral magnetic conductivity; the first model is Einstein grav- 
ity with U(1)l x U(1)r Maxwell theory, and our second set-up is based on the 
Sakai-Sugimoto model in a deconfined and chiral symmetry restored phase. While 
the former takes into account full back-reaction while the latter not, the common 
feature is an important role played by the appropriate 5-dimensional Chern-Simons 
term corresponding to the 4-dimensional axial anomaly. 
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1 Introduction 



For certain strongly coupled dynamics of gauge theories, gauge/gravity correspondence 
has become a useful method to study the problems, alternative to the conventional tech- 
niques such as perturbation theory. Although it was originally developed in the case of 
large color N c and strong t'Hooft coupling A = gy M N c limit, it has given us a consid- 
erable amount of new insights on generic strongly coupled gauge theories, and in many 
occasions its predictions capture, at least qualitatively, right physics for otherwise difficult 
non-perturbative phenomena. Of particular interests are of course possible applications 
to QCD. As QCD coupling runs to a large value at low energy, it is a logical hope that 
some low energy QCD phenomena which are hard to be explained by other means may 
have explanations in holographic QCD. It is at least worthwhile to study the problems 
in the framework to see what it pre /post- diets and also to compare with other known 
methods. 

Not only low energy QCD but also finite temperature deconfined quark-gluon plasma 
has been an active area of applications of gauge/gravity correspondence. The main moti- 
vation is the experimental finding at RHIC indicating a strongly coupled phase of quark- 
gluon plasma. Moreover, one may expect more sensible connections between finite tem- 
perature phases of different gauge theories and the QCD plasma, because at finite tem- 
perature fermions and scalar bosons get effective masses and become less relevant than 
the universal gauge field dynamics. Hydrodynamics would be a right place for searching 
for some universality as it is describing precisely the long wavelength transport dynam- 
ics for which these massive modes may decouple. There has been an enormous amount 
of recent works studying hydrodynamics and transport coefficients in the gauge/gravity 
correspondence ( see refs.[U[2] for reviews), although it still seems to remain as a fruitful 
area of further research. 

In this work, we will study one more example to the plethora of holographic QCD 
applications : the computations of chiral magnetic conductivity [31 H] at finite frequency. 
This is motivated by a recent work of Kharzeev and Warringa [I] which computes the 
time-dependent chiral magnetic conductivity in 1-loop perturbative QCD, aiming at a 
weakly coupled phase of QCD plasma at very high temperature. Chiral magnetic effect 
is the phenomenon where an electromagnetic current is induced parallel to the applied 
magnetic field in the presence of non-zero chiral density (i.e. unbalance between positive 
and negative helicity particles), and it is one kind of chiral-anomaly originated effects. 



1 



Chiral magnetic conductivity is the proportionality coefficient of the induced current to the 
magnetic field. In a QCD plasma such as the RHIC experiment, finite local chiral density 
can be generated by sphaleron fluctuations, and a large magnetic field may also appear 
in off-center collisions along the direction of angular momentum. Consequently, chiral 
magnetic effects may play some role in the subsequent plasma dynamics, and this was 
the motivation of the above authors. However, as RHIC plasma seems to become rather 
strongly coupled shortly after the collision, one needs other methods to complement the 
previous weak coupling calculation. Gauge/gravity correspondence would be a worthwhile 
try for this purpose. Lattice simulations for chiral magnetic effects for static but arbitrary 
magnitudes of magnetic field are given in refs.[5]. 

In fact, the possibility of using holography for computing chiral magnetic conductivity 
was first pointed out by Rebhan, Schmitt and Strieker in ref. [6]. Moreover, Son and 
Surowka recently computed a similar quantity (they call in the static magnetic field 
case [7]. See also the work by Lifschytz and Lippert [8] for other interesting phenomena 
related to chiral anomaly in the holographic set-up. However, it seems that there has 
been no study on the frequency dependent behavior of the chiral magnetic conductivity, 
which would be relevant in hydrodynamic simulations of the RHIC experiment. This will 
be our main focus and results in this paper. 

As the precise holographic model of large N c QCD is not known up to now, our objec- 
tive is to set-up a couple of consistent holographic frameworks to compute time-dependent 
chiral magnetic conductivity, and to provide the results based on those. Our first holo- 
graphic model is the 5-dimensional Einstein gravity with a negative cosmological constant 
coupled to two U(l) Maxwell fields; U(1)lX U(1)r. As should be clear from the notation, 
these two 5-dimensional gauge theories correspond to the 4-dimensional U(1)l x U{1)r 
chiral symmetry that we are interested in. We will consider an exact Reisner-Nordstrom 
black-hole solution charged only under the axial U{\)a = U(l)i — U(1)r to represent a 
finite temperature phase with non-zero chiral densitjo The " electromagnetism" , of which 
we will turn on an external magnetic field and also read off the induced current, is the 
vector part U(1)em = U(1)l + U(1)r. The second holographic set-up that we will study 
is based on the Sakai-Sugimoto model [10], which seems to be closer to the realistic QCD 
in quenched approximation. We will consider its deconfined, chiral symmetry restored 
phase with one flavor Np = 1, whose chiral symmetry is also U(1)l x U(1)r. The effective 
5-dimensional U(1)l x U(1)r gauge theory is a Dirac-Born-Infeld action and the back- 



1 See refs.[9] for previous studies of using this solution. 
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reaction of it to the background geometry is consistently neglected in quenched/probe 
approximation. This should be contrasted to the first holographic model where the full 
back-reaction is taken into account. 

In both set-ups, the common feature is the 5-dimensional Chern-Simons couplings for 
U(1)l x U(1)r gauge theory living in the holographic 5-dimensional bulk, corresponding 
to the 4-dimensional chiral anomaly. As this is more or less dictated uniquely by the 
anomaly structure, it is a universal feature, and the results directly related to it should 
be taken as robust. The only difference between the two models is the details of the 
background metric and how to treat the back- reaction of the finite axial/chiral density. 
One can take our set-ups as two prototypical examples whether or not we consider the 
back-reaction of the axial/chiral density. 

Note added : Shortly after this paper, there appeared an interesting observation in 
ref . [11] regarding the correct identification of holographic currents in the presence of 5D 
Chern-Simons terms, which is relevant in our computation of chiral magnetic conductivity. 
Ref. [1 lj computed zero frequency chiral magnetic conductivity in the Sakai-Sugimoto 
model taking into account additional contributions coming from these modifications. We 
will briefly summarize these modifications here, that will correct our currents we used in 
the text by a constant, frequency-independent shift. At the end, we will also mention a 
few puzzles that still remain even if we take this modification into account, which should 
be resolved in the near future. 

Firstly, as ref.pj] observed, one easily derives that the variation of 5D Chern-Simons 

term in the 5D action gives us additional contribution to the currents, 

]\r e ff at N e ff N 

A cs r L = -^-^e^(A L UF L ) pa , A cs r R = +^^e^(A R UF R ) pa , 

(1.1) 

which should be added to our currents in the text according to AdS / CFT dictionary. Note 
that the gauge fields appearing on the right-hand sides are external UV boundary fields, 
without any component of subleading dynamical piece. This is one characteristically 
different property of these additional contributions, which are of local, contact-term type, 
compared to the one from the subleading piece in the text. The total currents would be 
then the sum of the two if we accept these modifications. The above gives to the EM 
current Jem = e(jL + 3r) an extra piece 

AcsJem = * 2 V ^ CT ((A a )„(F EM ) pa + (A EM ) v (F a ) pa ) , (1.2) 
where the external EM and axial gauge potentials are defined by Al = eA E M — A a and 
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*C8*=-1SWN.) • (1.3) 



Ar = cAem + A a . In our situation, we turn on the time component of axial potential 
(A a )o = fj, a as a chemical potential for the axial charge, and also an external EM magnetic 
field B EM = (F EM )i2 along x 3 direction, to define chiral magnetic conductivity. From the 
above, the modification of chiral magnetic conductivity would therefore be 

Note that due to the local nature of the above modifications, the shift is simply a constant 
without any momentum or frequency dependence of the probe field Bem- Therefore, our 
plots can simply be shifted by this constant amount without a need for re-computations. 

This is one story, while ref. [UJ went further to propose an interesting observation. 
They realized that even after the above modification, the EM current is not strictly 
conserved in the presence of external axial potential A a , that is, one can find from the 5D 
equations of motion that 

N eff N 

^^V^™^ • (L4) 
The details can be found in ref.[UJ, but we only mention that they added an additional 
local counter-term, called Bardeen-term, in the regularized holographic effective action, to 
remedy this non-conservation. This additional counter-term, which can be added on the 
UV boundary as a different holographic renormalization prescription, gives us additional 
contribution to the current of similar type as above from the 5D Chern-Simons term. 
As can be easily expected, this contribution is also of local type and its contribution 
to the chiral magnetic conductivity is again a simple constant without any frequency 
dependence. 

If we choose to include this too, in total the induced EM vector current Jem = e(jL + 
jit) receives additional contribution 

p 2 Are// 7y 

AA/= ^ 2 C e^(A A UF EM ) P , , (1-5) 
which gives us a constant shift in the chiral magnetic conductivity by 

^ (Kreff; 



Act = [N F »N C ) , (1.6) 

that is precisely minus of the zero frequency value cr(0) we obtain in this paper, so that 
the zero frequency chiral magnetic conductivity in their prescription vanishes. As they 
pointed out, there is no a priori reason for vanishing chiral magnetic conductivity when 
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the vector current is strictly conserved. In fact, because the additional contribution is a 
frequency-independent shift, the real part of our resulting Figures I2|3ll5l and [6] are simply 
shifted down by Act, and there is non-zero chiral magnetic conductivity at finite frequency 
even in their prescription. As the magnetic fields relevant in RHIC experiments are time- 
dependent, chiral magnetic conductivity will still be at work, but with a different detailed 
prediction. 

Finally, let us mention a few seemingly puzzling aspects of the currents even after we 
take into account Chern-Simons contributions. First, combined with our resulting plots, 
the chiral magnetic conductivity goes to a constant value when u — > oo. Dynamically 
this doesn't make sense, because the medium cannot respond to the perturbation which 
is arbitrary fast. If we use only sub leading piece as we did in the text, we do get a nice 
damping when uo — > oo. This seems to indicate there might be something we are still 
missing at the moment. Another point, which is probably related to the first point, is 
that when we take a variation of 5D Chern-Simons action, one also gets a contribution 
from the IR boundary, which in our case is the horizon. It could be that one has to 
consider also the IR boundary contributions to resolve the first puzzle, which we hope to 
clarify in the near future. 

In summary, there indeed is an issue regarding what is the correct holographic current 
in the presence of 5D Chern-Simons term as we briefly reviewed current proposals. At 
the moment, arguably there seems no definite answer for that. However, the purpose 
in the present paper is to study frequency dependence of chiral magnetic conductivity, 
and because the present differences between different proposals are all constant shifts, our 
main results can easily accommodate the future resolving the issue. 

2 A quick review on physics of chiral magnetic effect 

A 4-dimensional field theory at finite temperature that gives rise to the chiral magnetic 
effect, such as chiral symmetry restored phase of QCD plasma, has the following basic 
ingredients 

• There are two chiral U(l) symmetries U(1)l and U(1)r, each having a non-zero 
triangle anomaly with the same magnitude but with opposite sign. Equivalently, if 
one weakly gauges these symmetries by coupling to non-dynamical gauge fields Al 
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and A R respectively, their current conservation laws are violated by 

N eff N 

d,f L = ^-^e^{F L ), v {F L U , 
N eff N 

d,f R = --^fe^{F R )» v {F R U , (2.7) 

where N F ^ is the effective number of flavors counted as fundamental representations 
of the color SU(N C ), and Fl,r are field strengths of Al,r. 

• Turn on a finite chemical potential for the axial U(1)a whose current is 

3a = -3l + 3r , (2.8) 

while keeping the system neutral under the vector "electromagnetic" U(1)em 

Jem = e(j L +j R ) , (2.9) 

where e is the electromagnetic coupling constant. In the language of weakly-gauging 
symmetries, these correspond to 

A L = eA EM -A A , A R = eA EM + A A . (2.10) 

Note that U(l) em is anomaly-free and it is consistent to let Aem be a real dynamical 
gauge theory as the Nature does. 

• Apply a homogeneous but possibly time-dependent magnetic field of Aem along say 
x 3 direction, 

B% M = (F EM )i2 = S(w)e- <w * , (2.11) 

which should be treated as an external perturbation to the system. Then the chiral 
magnetic effect induces the electromagnetic current j EM parallel to B EM 

3%m = 3^)^ = cr(u)B(u)e-^ . (2.12) 
The <j(oj) is the chiral magnetic conductivity [21 H]. 

An intuitive explanation on the microscopic origin of this phenomenon was given in 
ref.[l]. For simplicity, consider a free massless one flavor of quarks {qLiQu) with unit 
electromagnetic charge. We emphasize that masslessness is important to have chiral 
symmetry. Upon quantizing qi, one gets a particle of negative helicity (meaning that its 
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spin is opposite to its momentum) as well as its anti-particle of positive helicity with a 
negative electromagnetic charge, and vice versa for q^. Due to the Wigner-Eckart theorem, 
the magnetic moment of an elementary particle should be proportional to its spin, and for 
positively charged particles/anti-particles the magnetic moment is in the same direction of 
the spin, while for negatively charged ones it is opposite to the spin. Now imagine applying 
an external electromagnetic magnetic field, then the magnetic moments of particles/anti- 
particles will tend to align along the direction of the applied magnetic field. Because 
the magnetic moment, the spin, and the momentum are correlated with each other as 
described above, one can easily deduce the following pattern of responses 

• Positively charged particles as well as negatively charged anti-particles from q L tend 
to move in reverse direction to the magnetic field. Let's denote them as (q\/ 2 , 9+1/2) 
where the upper index represents the charge and the lower the helicity. 

• Positively charged particles as well as negatively charged anti-particles from q^ 
tend to move towards the same direction as the magnetic field. We denote them as 

(9+1/2 > 9-1/2)- 

Then, having a finite chemical potential for the axial U(1)a symmetry which might be 
achieved by local sphalerons means that the number of positive helicity states is larger 
than the number of negative helicity states 

N (9+V2) + N « 1/2 ) > N (qt l/2 ) + N (qZ l/2 ) . (2.13) 

In conjunction with the above discussion, observe that the left-hand side of the above 
inequality induces a positive electromagnetic current along the magnetic field, while the 
right-hand side would contribute to a current in opposite direction to the magnetic field, 
so that the above inequality tells us there would be a net positive electromagnetic current 
induced along the magnetic field : this is the chiral magnetic effect. 

Although it is not absolutely necessary, the electromagnetic neutrality that we require 
for simplicity implies that 

N (£1/2) " N (£1/2) =~( N (€1/2) ~ N (9:1/2)) , (2.14) 

where the left-hand side is proportional to the chemical potential for qL (or U(1)l) and 
the right-hand side is the minus of the chemical potential for q^ (or U(1)r), which implies 
jjji = = jtt a with ha being the axial chemical potential we turn on in the background. 
It might be an interesting future direction to generalize our computations to the cases 
with non-zero electromagnetic charge density too. 
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3 Retarded response (Green's) function in Eddington- 
Finkelstein coordinate 



Before presenting holographic calculations of chiral magnetic conductivity, let us make a 
short digression to explain our method of obtaining retarded response function in a black- 
hole background. There is by now a well-established procedure to compute any retarded 
response function in gauge/gravity correspondence in linear response approximation [T2] , 
and our results will also belong to this category. As we are interested in the finite frequency 
u, the method based on derivative expansions such as those of ref.|13j. although it is fully 
non-linear, is not suitable for our purpose, as one necessarily truncates the expansion 
at some finite order of derivatives and the results lose its validity at high u. On the 
other hand, the linear response approach can be used at any frequency uj while it requires 
the driving source to be small for the linear approximation to be valid. It seems both 
approaches have pros and cons. 

In linear response approach, one perturbs the system by an external source B coupled 
to an operator J of the theory, and the first-order perturbation theory tells us that the 
response in the expectation value of the operator (J) is given by the convolution of the 
source with the retarded Green's function 



Typical interests, including our present work, are therefore the computations of the re- 
tarded response function Gr. In the language of quantum mechanics, the retarded re- 
sponse function Gr has an expression 



and it is one kind of 2-point Green's function in Minkowski signature. Because gauge/gravity 
correspondence is typically formulated as a prescription to compute precisely the Green's 
functions or correlation functions, one might hope to compute Gr rather easily by simply 
applying the suitable gauge/gravity dictionary. However, this "suitable" dictionary turns 
out to be rather non-trivial in Minkowski signature, especially retarded Green's function 
in the presence of black-hole horizon, issues being whether or not one should treat the 
horizon as boundary, etc. The issues have been settled by now, and one has a definite 
well-defined way of computing Gr holographically in linear response theory |12j . 

However, looking back the original motivation of studying the causal response in the 
presence of driving external source, the quantum mechanical expression (I3.16P which is a 




(3.15) 




(3.16) 
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language of 4- dimensional field theory side in fact seems unnecessary. What one is inter- 
ested is simply the resulting ( J(x)) causally responding to the external source B(x) in the 
Minkowski evolution of the 5-dimensional holographic dual theory. The fact that Gr has 
a field theory interpretation of 2-point correlation function is not needed to find the an- 
swer, if one can directly solve the Minkowski dynamics in 5-dimensions in the presence of 
external source B(x), because the expectation value (J{x)) is also encoded in the resulting 
5-dimensional solution as a normalizable mode in the near boundary expansion according 
to the holographic renormalization [H]. In other words, once we accept the results of stan- 
dard holographic renormalization a la ref. [14J, which states that the expectation values 
can simply be read off from the near boundary expansion of the 5-dimensional fields even 
in the presence of black-hole horizon and in Minkowski signature, we can by-pass the is- 
sue of calculating retarded Green's function, by simply solving the Minkowski equation of 
motions with "physically obvious" in-coming boundary condition at the black- hole hori- 
zon. With the UV boundary condition also fixed by the given source B(x), this uniquely 
determines the 5-dimensional solution and from it one can directly read off (J(x)). The 
important question is whether this way gives us the same results consistently to those 
obtained by the well-established method of computing retarded Green's functions. For a 
massless scalar, it is checked to be true, and we conjecture it is always true. 

Turning to a practical side, we will adopt the Eddington-Finkelstein coordinate system 
for a given background black-hole geometry in which the metric takes a form 

3 2 

ds 2 = -r 2 V(r)dt 2 + 2drdt + r 2 £ (dx*) , (3.17) 

i=l 

where the horizon is located at the position V{ru) = 0. The motivation is that the 
in-coming boundary condition at the horizon can be easily fulfilled in this coordinate by 
simply imposing only regularity at the horizon on the solutions. To see this more clearly, 
consider a mode with frequency to along the time-like Killing vector & 

<Kr,t)~e**f(r) . (3.18) 

Near the horizon where V(r) ~ 0, the metric becomes 

ds 2 ~ Idrdt + ds 2 T , (3.19) 

where ds\ is along the transverse directions x l that is not important in the discussion. 
Therefore, the (r, t) coordinates become two local null directions near the horizon, and in 
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Figure 1: A schematic description of the horizon in the Eddington-Finkelstein coordinate. 

terms of the more standard local flat coordinates, say (X ,^ 1 ), they are written as 

t = -L(x» + X>) , T + , (3.20) 

so that ds 2 ~ — (dX ) 2 + (dX 1 ) 2 + ds^. In writing the above, we re-parameterize r such 
that r = is the location of horizon. See the Figure [I] for a schematic definition of 
(X ,^ 1 ), and the region of X° > X 1 is inside the horizon while X° < X 1 corresponds to 
the region outside the horizon. Hence, the mode (13.1 8ft on which we impose only regularity 
at the horizon r = r# has a near horizon behavior 

<P ~ e-^f(r H ) = e- % ^ X ° +Xl )f{r H ) , (3.21) 

which is automatically in-coming towards the region inside the horizon. Note that the 
conclusion doesn't depend on the signature of u. 

Another way of understanding this is to go back to the more conventional form of the 
black-hole metric 

ds 2 = -r 2 V(r)dtl + ^- )+ d4 , (3.22) 

achieved by a relation 

rr Apr 1 

U=t -Lww) ■ (3 - 23) 

Note that the time-like Killing vector retains the same expression 

1 = 1: • (X24) 
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and a generic mode with frequency u would be 

0(r, Q = e~ iuJt * U (r) = e~ iuJt (e™ & M /» (r)) = e"^7(r) , (3.25) 
so that by definition 

f m [r) = e~ luJj ™ (r'pv(r') /( r ) . (3.26) 

Near the horizon the phase factor in (13.261) is precisely what one would need to make 
(j)(r, £*) to be an in-coming wave, so that by simply requiring only regularity of f(r) 
at the horizon the mode becomes in-coming automatically ; in other words, working in 
the Eddington-Finkelstein coordinate is equivalent to giving a preference to the in-coming 
modes, and the out-going modes in (r, £*) coordinate would look singular ia the Eddington- 
Finkelstein coordinate. Reversing the logic, an inverse Eddington-Finkelstein coordinate 
where 

ds 2 = -r 2 V(r)dt 2 ^ - 2drdU* + ds 2 T . (3.27) 

should be useful for picking-up only out-going modes and for computing advanced response 
(Green's) functions. 



4 Holographic model I : Einstein plus U(l) 2 

The simplest holographic model that implements the symmetry structure of the section 
2 would be a 5-dimensional Einstein gravity with negative cosmological constant plus a 
U(1)l x U(1)r gauge theory corresponding to the chiral symmetry in 4-dimensions. In 
addition, the 4-dimensional triangle anomaly of U(1)l X U(1)r will be reflected as a 5- 
dimensional Chern-Simons term in the holographic model. The minimal Lagrangian is 
then 

(16nG 5 )£ = R+12-±(F L ) MN (F L ) MN - \(F R ) MN (F R ) MN (4.28) 
+ i ^ e A ™«^(A i ) M (F L )^(F L )Q fi -(^) M (F fl ) W p(F i? ) Q ^ , 

where we normalized the cosmological constant to have a unit radius for AdS§ for simplic- 
ity. The capital letters M, N, . . . denote 5-dimensional indices while Greek letters /i, u, . . . 
will be reserved for 4-dimensional coordinates. The epsilon symbol in the above is purely 
numerical as we factor out yf—gi explicitly, and our convention is e rfiual3 = ^ va ^ . 
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It is not difficult to relate the value of k with the 4-dimensional anomaly coefficient 
[?]. The equation of motion for Al is 

Vn{F l ) mn - — ^=== e MNPQR (F L ) NP (F L )Q R = , (4.29) 

and in the pure AdS§ vacuum in Poincare coordinate 

ds 2 = -r 2 dtl + ^ + r 2 j2(dx 1 ) 2 , (4.30) 
r i=i 

the M = r component becomes 

d lt (F L )^-^^(F L ) lv (F L ) aP = . (4.31) 

According to AdS/CFT dictionary, the 4-dimensional U{1)l current that A L corresponds 
to is given by 

•^-lek^'' 3 ™'" ' <4 ' 32) 
while the external weakly gauging potential that couples to j£ is the boundary value of 
A L {r = oo), so that by taking r — > oo limit of (I4.3ip . we have 

d^ L = -^-/ uaP {F L ), v {F L U , (4.33) 

where Fl now represents the 4-dimensional weakly gauging potential coupled to ji- In 
comparison to (12. 7p . we thus have 

K = - 2 -^(N^N c ) . (4.34) 

Because we are going to turn on the axial chemical potential, it is more convenient to 
work in terms of A em and A a related to Al,r by ( 12.1 Oft . 

Al — &A em — A a , A R = eA em + A a , (4.35) 

where these fields now represent 5-dimensional gauge fields accordingly. The action den- 
sity then takes a form 

(16ttG 5 )£ = R+l2-^(F em ) MN (F em ) MN - ]-(F a ) MN (F a ) MN (4.36) 
€ mnpqr f 3e 2 ( A g ) ^ (F em ) NP (F em )Q R + (A a ) M (F a ) NP (F a )Q R ) , 
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whose equations of motion are 

Rmn + [ 4 + — {F em ) 2 + - (F a ) 2 ] g MN — e 2 (F em ) PM (F em ) P N — (F a ) PM (F a ) P , 



6\ cu t / • r 1 
o 



^N{F a ) MN + . , e MNp QR (e 2 (F em ) NP (F em )Q R + (F a ) N p(F a )n R ) — , 

V N {F em ) MN + .r^e^^iFJMFemhn = . (4.37) 

The normalization of the A em kinetic term looks a little unconventional, but in this way 
the UV boundary value of A em couples to the EM current j em = e{ji + jn) with unit 
strength, and it is more convenient in this sense. 

There is an exact AdS Reisner- Nordstrom type black-hole solution to the above equa- 
tions of motion with only an axial charge being turned on, and we will use this space-time 
as a background representing a finite temperature plasma with an axial/chiral chemical 
potential. The explicit form of the solution in Eddington-Finkelstein coordinate is 

3 

ds 2 = -r 2 V(r)dt 2 + 2drdt + r 2 Y / (dx 1 ) 2 , 
A a = f — g — j dt , A em = , 



H 



nr) = 1-jS + f? • (4-38) 

where the location of horizon is the largest root of V{th) = 0, and the axial/chiral 
chemical potential can be identified as 

(4-39) 

Our task is to perturb the above solution by a (small) oscillatory external EM magnetic 
field Bl m along x 3 , and solve the resulting linearized equations of motion of (14.371) to read 
off an induced EM current j 3 m along x 3 direction. Holographic renormalization tells us 
that jf m is encoded in the near UV boundary expansion of the resulting solution of A em 
by 

e 2 

3em = TVT limr 2 (A ero ) 3 ■ ( 4 - 40 ) 

47T(jr5 r ^°° 

Note that the factor e 2 is from our convention of the 5D kinetic term for A em in the action 
density. 

In this linearized order, it is easy to see that the gravity fluctuations 5 gMN and the 
fluctuation of the axial gauge field 5A a in fact do not couple to the electromagnetic 
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fluctuation d~A em , so that the problem is much easier than naively expected. For example, 
in the second equation of (14.37ft it is consistent to turn off 5A a = as 5A em 's contribution 
would be only quadratic. The same is true for 5gMN i n the Einstein equation, so that for 
our purpose we can simply let Sqmn = 5A a = and just consider the equation 

V N {F em ) MN + -^e MNp Q R {F a ) NP {F em)QR = , (4.41) 

with A a being replaced by the background value in the above solution. This decoupling 
between SA em and 8gMN has a physical interpretation : to linear order, external electro- 
magnetic perturbations do not introduce additional energy-momentum, especially there 
would be no momentum flow induced at first order. For example, in the microscopic 
picture for chiral magnetic effects in section 2, if we keep electrical neutrality by 

N {V-1/2) ~ N (q+ 1/2 ) = ~ (N (ql 1/2 ) - N ( qZl/2 )) , (4.42) 

the left-hand side is proportional to the chemical potential /jLl for gx, and vice versa for 
the right-hand-side and qR. Given a temperature, the thermal physics of qL would be 
precisely mirror to the physics of qR, so that we can expect that the number of iV (q^ 1 ^^j 
from q L is precisely equal to iV (qZi/ 2 ) from q R in the situation of hl = —[Ir. Similarly, 
N (q+1/2) — N (qXi/2) ■ This gives us at the end 

N (qt 1/2 ) + N (q- 1/2 ) = N (g+ /2 ) + N (qZ 1/2 ) . (4.43) 

However, observe that the right-hand side is the total number of particles/anti-particles 
that move along the magnetic field, while the left-hand side is the total number that move 
opposite to it, so that we conclude there would be no net momentum flow induced along 
the magnetic field. This is the physics reason behind the decoupling between SgMN and 
SA em , which should be universal for any EM perturbations in neutral systems. 

This feature would be absent in the set-up having only one U(l) symmetry with a 
triangle anomaly, as recently studied in ref.[7]. In this case, we have, say, q^ only without 
q R . Turning on a chemical potential for qL, one would have an excess of q\/ 2 over qZi/ 2 
; N(qZ 1 / 2 ) > N(qZ l i 2 ), and an external magnetic field would move them together in the 
opposite direction to it, and would subsequently induce a net negative chiral magnetic 
current. This effect corresponds to £b in ref.[7] for a static external magnetic field. How- 
ever, as the particles/anti-particles from qL move together in the same direction, without 
having a compensating flow by qR, there will be a net momentum flow introduced by the 
effect. This complicates choosing a right Landau frame on which there shouldn't be any 
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momentum flow by definition. This might explain the complicated second term of in 
ref.0. 

Back to our main purpose, the linearized equations for A em in components are 

M = r : V N F rN = , 
M = t : V N F tN = , 

M = i : V N F iN + ^-e ijk F jk = , (4.44) 

where we omitted and will omit 5 and subscript em in 5A em for clarity in the following. 
The non-vanishing Christofel symbols are 

Kt = n r = -T\ t = - l -d r (r 2 V{r)) , VI = \r*V{r)d r (rV(r)) , 

r;,=r} r = ^ , I-=-rV(r)^ , r^ = -n% , (4.45) 

from which one arrives at the following explicit equations 

(d t F tr ) + 1 + V(r) («9,F„) = , 

d r (r 3 F tr ) - r (^F ri ) = , (4.46) 
8 r (rF tl + rV(r)F„) + r (<9 t F„.) + 1 (6^) - ^e^ fc F ifc = . 

We should fix our solution ansatz in accord to our purpose of computing time-dependent 
chiral magnetic conductivity. First, one can always work in the radial gauge A r = 0. 
Moreover, as we wouldn't expect any net EM charge density appears in response to the 
homogeneous magnetic field, one can also adopt the ansatz A t = 0. We will see the con- 
sistency of this ansatz later as it will consistently solve the full equations of motion. We 
assume a definite frequency u> for every field 

A i = A i {r,x i )e- iut , (4.47) 

and we require homogeneity along the x 3 direction in which we turn on the magnetic field 
B 3 = Fi2, so that we simply drop 83 in the equations of motion. Again, the consistency 
of these can be checked by the full equations of motion that our solutions will solve. As 
F tr = in our ansatz, the first two equations in (14.461) are uniquely solved by 

ftAi + d 2 A 2 = . (4.48) 

At least we should have a non-zero F u = d\A 2 — d 2 Ai, so that we need to give a finite 
wave-number k± = (ki,k 2 ) to the transverse x — (x - X ) coordinate, although we will 
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take a homogeneous limit k± — > at the end to get a finite value of chiral magnetic 
conductivity with only frequency dependence. These steps finally give us the following 
ansatz 

A x = k 2 f{r)e- iuJt+lk± - x± - , 
A 2 = -k 1 f(r)e~ iu ' t+ik± - x± , 

A 3 = g(r)e~ iwt+ik ^ x± . (4.49) 

The remaining equation to solve, i.e. the third equation in (14 .46 p . gives us the following 
set of two equations for f(r) and g(r), 

d r (-iuirf + r 3 V{r) (d r fj) - lur (d r f) - -k 2 J ' ' 2h ' Q 







d r (-iuorg + r 3 V{r) (d r g)) - iuor (d r g) - -k\g + % -^-k 2 J = , (4.50) 



1,, H2kQ . 
~ k ±9 + — 3 

With these being solved, the full equations of motion are satisfied within our ansatz. 

It is important to determine the right boundary conditions for the radial profiles 
f(r) and g(r). As discussed in section 3, simple regularity at the horizon r = r^ will 
be enough to implement the in-coming boundary condition that is suitable for a causal 
retarded response. At the UV boundary r — ► oo, we need to have an external magnetic 
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field B 3 = F\n. Because of 
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-ik 2 ± f(r)e- iujt+ik± - x , (4.51) 



this dictates that /(oo) 7^ 0. On the other hand, we shouldn't get any external perturba- 
tion from A3, which imposes the normalizable boundary condition on g(r), 

9( r )^°(^) » r^oo . (4.52) 
In fact, the induced EM current along x 3 will be obtained by 

jL = ^ lim^e'^^ , (4.53) 
so that our desired chiral magnetic conductivity will be computed as 

a(u,k ± ) = -^ r \im^f\ . (4.54) 
v ' 4irG 5 r-oo k 2 ± f(r) v ' 

One can easily convince oneself that the equations (I4.50p with the above boundary condi- 
tions pose a well-defined procedure of calculating a(u>, k±). One should resort to numerical 
jobs to proceed further, however. 
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Focusing on a special case of homogeneity limit k± — > 0, one indeed gets a finite value 
of chiral magnetic conductivity. Expanding f(r) in power series of k±, 

f = fo + k\h H , (4.55) 

one finds that g(r) should start its expansion from 0(k±) for consistency, 

g(r) = k 2 ±9l + ■■■ . (4.56) 

This is because the homogeneous differential equation 

d r (-iujr ■ +r 3 V(r) (0 r -)) - iojr (d r -) = , (4.57) 

that go would have to satisfy has a unique solution regular at r = r# up to an overall 
multiplication factor, and its UV boundary value at r — > oo is in general finite unless one 
encounters a quasi-normal mode, but this happens only with a complex oj with a negative 
imaginary part. However, this would contradict to our boundary condition (14.521) . and 
one necessarily has g = 0. 

Inserting the above series expansions into (14.501) . one gets 



d r (-iurfo + r 3 V(r)<9 r /o) - iwr (d r f ) = 

H2kQ 



d r (-iur 9l + r 3 V(r)d r9l ) - tur (d r9l ) + = , (4.58) 

which is enough to obtain the homogeneous limit of chiral magnetic conductivity 



<7(w) = <7(w, k x ^ 0) = ^— lim . (4.59) 

Numerical job for computing ct(c<j) with (I4.58P is much simpler than the general case of 
o~(u, k±) : One can first solve fo from the first equation, and then use this to obtain g\ 
subsequently. See Figure [2] and Figure [3] for the numerical results of cr(u) for several differ- 
ent axial chemical potentials. In the Appendix, we outline an easy method of producing 
our numerical results. 

Zero frequency limit uj — ► 

It is useful to check our computation for the static limit u — > as it has been argued 
that the static limit is universally dictated by anomaly irrespective to the details of the 
dynamics. The equation for / simplifies to 

d r (r 3 V(r)d r f ) = , (4.60) 



17 




Figure 2: Frequency dependent chiral magnetic conductivity cr(u) for various axial chemical 
potentials. The solid line is the real part of cr(u) while the dashed one is the imaginary part. 



whose unique regular solution at r = r# up to normalization is a constant : f = 1. The 
subsequent equation for gi 

d r [r 3 V(r)d r gi) + —^- = , (4.61) 

with the regularity at r = r# as well as the normalizability C I4.52I) at UV determines the 
unique solution 

* w = (TO fe-i) • <« 2 > 

Using the explicit form of V(r), we in fact only need the fact that V{r) — > 1 as r — > oo 
to find the near boundary behavior of g± as 

9i{r) -> ^5^-^ H , r^oo , (4.63) 

so that the chiral magnetic conductivity from (I4.59P at zero frequency is 

^ - o) ^ (^) ■ <« 4 > 

where we used (I4.34p 

K = -??*( N *ff Nc ) , (4.65) 

37T V ' 

as well as the expression for the axial/chiral chemical potential ( 14.391) 

A«a = %■ , (4.66) 



ju/T=l /i/T=10 




Figure 3: More results for other values of axial chemical potentials. 



in the last equality. This is indeed the right answer. 

5 Holographic model II : The model of Sakai and 
Sugimoto 

Our second model for a holographic calculation of time-dependent chiral magnetic conduc- 
tivity is the deconfined and chiral symmetry restored phase of the Sakai-Sugimoto model 
at finite temperature. This model is supposed to describe, at least qualitatively, the large 
N c limit of QCD with a small number of fundamental quarks in quenched approximation, 
and in this sense would be more realistic than the somewhat arbitrary set-up in the previ- 
ous section. For simplicity, we will focus on a single flavor Np = 1 case, which means that 
we have a single D8 and D8 probe brane pair embedded into a known black-hole solution 
for a deconfined phase [15]. The chiral symmetry U(1)l and U{1)r that we are interested 
in live on the world- volumes of these probe D8 and D8 branes respectively. For the chiral 
symmetry to be restored/unbroken, we confine ourselves to the phase where each of these 
branes meet the black-hole horizon and do not meet with each other. See the Figure H] for 
a schematic picture. This also implies that the leading order dynamics on each 8-branes 
are independent of each other : assigning U{1)l to D8 and vice versa for U{1)r and D8, 
the total action will simply be a sum of the two 8-brane world volume actions 

Stat = S m {A L ) + S m {A R ) , (5.67) 
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X 4 



Black— hole horizon 

Figure 4: A schematic picture of the Sakai-Sugimoto model in its deconfined and chiral sym- 
metry restored phase. 



where 

Sm/BS = -I* J ^ \/ det (9* + 2 ^ F ) T ^ ( / F* R AAAFAF , (5.68) 

with n p = (27r)~ p lj (:P+1 ' . Note that the Chern-Simons coupling for D8 has the opposite 
sign to that of D8, and we show only the relevant term with the F RR 4-form Ramond- 
Ramond field strength as our background carries only this Ramond-Ramond field. 

The 10- dimensional black-hole solution that provides a background corresponding to 
a deconfined phase in Eddington-Finkelstein coordinate is [15] 



ds 2 = (^J y-f(U)de + (dx l ) 2 + -^— Yjdxlj +2dUdt+ (-) [7 2 dft 2 
F 4 - = ^#^e 4 , e^^V , ^Vol(^) 
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V4 " V/?, 

i? 3 = ng,N c l 3 a , f(U) = l-(jf) , (5.69) 

where the temperatute j3 = ^ is related to the location of the black-hole horizon Ut by 

47T / i? 3 \ ^ 

"=-ry • (5 ' 70) 

and has a period (27rZ s ) which is not essential for our purpose though. The F RR is 
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normalized in such a way that 

F* R = N C , (5.71) 



(2^/ s ) 3 Js* 

and this convention conforms to those in our writing of 8-brane actions (15.681) . Each 
8-branes of our interest spans the coordinates (U, t, x\ Q^) of total 9-dimensions sitting at 
a point in £ 4 , but dual QCD dynamics should be insensitive to the modes inside S 4 , so 
that we will assume homogeneity along Q4 directions and integrate over them from the 
very first stage. After that, the resulting 5-dimensional action takes a form with a minor 
computation 

Sds/ds = -CBJ J dx A dUU^det{g* + 2irP s F) T^J dx A dU e MNPQR A M F NP F QR , 

(5.72) 

with a definition of 

C = 

3-2 5 -7riT£fz; 



N2 

C = —r^r , (5.73) 



for later convenience. 

Our objective is to first construct a background solution having a finite chemical 
potential fi a for the axial current 

A a = l -{-A L + A L ) , (5.74) 

while keeping the system neutral under electromagnetism for simplicity 

eA em = i (A L + A R ) , (5.75) 

and then to consider a small external perturbation of electromagnetic magnetic field to 
see its retarded response in the current. This means we turn on a chemical potential 
fj>L = —fJ-a for the U{1)l symmetry on the DS brane, and similarly = /z a for U(1)r on 
the D8 brane. See ref.[T6l [TTJ for turning on iso-spin chemical potentials in the model. 
To find the corresponding background solution, it is enough to keep F tu component only 
in the action of the 8-branes, and a simple calculation gives 

SmfM = ~cj d'xdUUi^l - (2nl^(F tu y , (5.76) 

whose solution is easily integrable to be 

F tu = , ±a ^= , (5.77) 
iJU 5 + (2< 2 ) 2 a 2 
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with an integration constant a > 0, where the upper sign is for Al on the D8 and vice 
versa for Ar on the D8. One can always work in the radial gauge A v = in which the 
above solutions can be further integrated for A t as 

MU) = =F f dU> - Q = , (5.78) 

^(rj') 5 + (27r/2) 2 a 2 

where we have imposed the condition A t (C7 = Ut) = at the horizon to fix an integration 
constant. The chemical potentials of U(1)l,r are then read off as the UV boundary values 
of the above A t respectively which relates a with fi a as 

f 00 2a p / 3 1 13 (2^ 2 ) 2 q 2 \ 

H a = / rfZ/ . = = r 2-fi — , -, — , 775 • (5.79) 

J Ut ^{U'f + (27r/ 2 ) 2 a 2 3C/| V 10 2 10 U t J 

We should remember this relation to identify a physical meaning of a. 

Having obtained the background solution, one next needs to find a linearized equation 
of motion for the vector part electromagnetism treated as a perturbation to the back- 
ground solution. One way of performing the analysis is to expand the 8-brane actions up 
to second order in fluctuations so that the linearized equation of motion can be derived 
directly from it. Using the series expansion 

y/det (1 + 5 A) = 1 + ^tr (6 A) + ^ [tr (5 A)] 2 - - [tr ((SA) 2 )] + O ((5A) 3 ) , (5.80) 

a straightforward but substantial amount of algebra finally produces the result 



S nl/m = cjd^xdu 



L{U) + \A{U) (SF tu ) 2 - B{U) (5F ti ) (5F m ) - l -C{U) (5F m ) 2 



- 1 -D(U)(5F t] )(5F tJ ) 



2 v ,v v , v v 2 



8tt 2 

where here A t (U) in the last term is 



d 4 xdU e^ k A t (U) (6F m ) (SF jk ) , (5.81) 



A t (U)=+f dU' -, a = , (5.82) 
J Ut J(U'f + (27r/ 2 ) 2 a 2 



for both D8 and D8 branes, because the sign difference from the Chern-Simons terms is 
compensated by having the opposite background chemical potentials, so that the fluctu- 
ation actions for the 8-branes are now identical to each other. The appropriate functions 
that appear in the above are 

L(U) = -Ul (l- (2nl 2 s ) 2 {F tu ) 2 f = -U 5 (U 5 + {2nl 2 s ) 2 a 2 y L2 , (5.83) 
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A(U) = (27ra 2 [/f (l-(27ra 2 (^) 2 )" f = (27ra 2 ^ 5 (^ 5 + (27rZ s 2 )V) f , 

3 3 

B(U) = (2^ s y (J) *ui (l - (2< 2 ) 2 (F tl/ ) 2 ) " 4 = (2vr/ 2 ) 2 (5) ' + (2< 2 ) 2 a 2 ) 1 

C(tf) = (27r/ 2 ) 2 /(f/)t/f (l-(2< 2 ) 2 (F w ) 2 ) 4 H2^) 2 /^(t/ 5 + (24 2 ) 2 «f , 
D(U) = (2VT/ 2 ) 2 (l - (2< 2 ) 2 (F tu f) 1 = (2< 2 ) 2 (J)' U 5 (U 5 + (2*# 'a 2 )"' 

For a qualitative understanding, the detailed form of the above functions are not essential 
except the fact that C(U) — > at the horizon U — > C/y, and their near-boundary behaviors 

A(17) - (2< 2 ) 2 ^l , B(U) -> (2< 2 ) 2 J Rl 1/ , 

C(17) -> (2tt/ 2 ) 2 f/i , D([/) -> (2nl 2 s ) 2 R 3 [H , U -> oo . (5.84) 

The structure of the above fluctuation action is in fact qualitatively identical to that in 
our first holographic model, which will be more manifest later. 

As the D8-brane action for fluctuations up to second order has an identical form to 
that of the Z)8-brane, the total action for the electromagnetic fluctuations would be simply 
a twice of (15.811) with the replacement 

SF M n — > e(5F em ) M N , (5.85) 

where we will omit the subscript em from now on for clarity, while we will leave 5. From 
this action with the near-boundary behaviors (I5.84p . a gauge/gravity dictionary for the 
current can be easily deduced as 

f em = 3e 2 C(27r/ 2 ) 2 Jim *7§ (<L4 3 ) , (5.86) 

which we will use later to find the induced electromagnetic current from the solution of 
5 A in response to an external EM magnetic field. 

It is straightforward to write down the equations of motion from (15.811) . 

d v (A(U)5F tu ) - B(U) (diSFm) = , (5.87) 
A(U) (d t 8F tu ) + B(U) (didFu) + C(U) (d t 5F m ) = , 

B(U) (d t 5F m ) + du (B(U)5F tl + C(U)8F m ) + D(U)d 3 (SFjA - -^-F tu e^ k 5F ]k = , 

whose structure is essentially similar to (14.461) in comparison as it should be expected. 
Subsequent analysis henceforth is quite close to the previous holographic model. Starting 
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from the ansatz 

SAy = k 2 f(U)e- iujt+ik± - x± , 

SA 2 = -kif(U)e~ iu;t+ik± - x± , 

5 As = g(U)e~ iuJt+ik ^ x± . (5.88) 

the full equations of motion are solved by the following two equations for f(U) and g(U), 

iN 

du {C{U) (d v f) - iuB(U)f) - iuB(U) (d v f) - D{U)k\f - ^F tu (U)g = , 

iN 

du (C(U) (dug) - iuB(U)g) - iuB(U) (dug) - D(U)k\g + ^F tu (U)k 2 ± f = , 

(5.89) 

with the boundary conditions that they are regular at the horizon U = Ut and g(U) 
should be normalizable at UV, 

g(U)~\ , U^oo . (5.90) 

U 2 

Using the current formula (15.861) . and the external magnetic field 

5F 12 = -ik 2 ± f(oc)e~ Mk± - x± , (5.91) 
the momentum dependent chiral magnetic conductivity will be computed as 

3 

a(^) = *3e 2 C(2< 2 ) 2 hm . (5.92) 

For the homogeneous limit k± — > 0, one expands 

f = fo + k 2 ± h + --- , g = k 2 ±9l + --- , (5.93) 

and /o and gi are solved by the equations 

d u (C(U)(d u f o )-iu;B(U)f )-iuB(U)(d u f Q ) = , (5.94) 

iN 

du (C(U) (dugi) - iujB(U)gi) - iuB(U) (d u9l ) + -^F tu (U)f = , 

where one first solves the former equation and then use that to solve the second equation. 
The homogeneous but frequency dependent chiral magnetic conductivity is then 

a(u) = i^C^f lim . (5.95) 
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Static limit uj — > 

We check the static limit of our formula for consistency. In oj — > limit, the equation 
for f is 

d u (C(U)d u f )=0 , (5.96) 

whose regular solution at the horizon is simply a constant, fo = 1, due to the fact 
C(Ut) = 0. The equation for g\ is then 

du(C(U)d ugi ) = -^F tu = du ^A t (U)) , (5.97) 
whose integration is 

ziV c A(£7) 

^ 1 = 4^^) ' (5 ' 98) 
where we have used the fact AtiUx) = so that the right-hand side is regular at U — Ut- 
Integrating once more, one gets 

where we already chose an integration constant to have a normalizable solution for g\. 
From the UV asymptotic behavior of C(U) 

C(U) -> (2tt/ s 2 ) 2 £/§ , U -f oo , (5.100) 

and also the previous identification of the chemical potential 

A t {oo)=fi a , (5.101) 

one finally arrives at 

^--M^m + - ' u — ■ (5 - 102) 

Inserting into (15.951) . one checks that 

a (uj = 0)=a = ^-N c , (5.103) 

which is the right result for a single flavor Np = 1 quark we are considering. 

To perform numeric jobs in the case of Sakai-Sugimoto model, one has to fix the 
parameters of the theory. First one can always work in the unit where 27r/ 2 = 1. From 
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Ai/T=0.01 p/r=o.i 




Figure 5: Time-dependent chiral magnetic conductivity <t{uj) for various axial chemical poten- 
tials in the Sakai-Sugimoto model with T = 200 MeV. The solid line is the real part of cr(u) 
while the dashed one is the imaginary part. 



the p meson mass and the pion decay constant, Sakai and Sugimoto determined the 
parameters [TOj [18] 

gl M N c ^ll , M KK = 0MGeV , (5.104) 

and we will take these values for illustrative purposes. For the temperature, we choose 
T = 200 MeV as a representative value relevant for the RHIC experiment^]. With these 
being fixed, there is no other ambiguity in the model. We present the numerical results 
of cr(u)) in Figure Eland Figure El 

6 Summary of results 

We set-up a couple of holographic frameworks for computing time-dependent chiral mag- 
netic conductivity. In the first model, we consider a full back-reacted Reisner-Nordstrom 
black-hole solution with only an axial chemical potential turned on, to study the induced 
electromagnetic current in response to a small, time-dependent magnetic field perturba- 
tion. Our second model is based on the more realistic model of Sakai and Sugimoto in 
its deconflned and chiral symmetry restored phase, but within a quenched/probe approx- 

2 According to ref.[15], the confinement/deconfinement transition happens at T c = M * K « 150 MeV 
in this model. See also ref.pjj] for further aspects. 
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Figure 6: More results in the Sakai-Sugimoto model with T = 200 MeV for other values of axial 
chemical potentials. 

imation. Both models give us qualitatively similar results for the frequency dependent 
chiral magnetic conductivity, which may be a useful complementary computation for the 
strongly coupled regime to the existing recent weak-coupling computation in perturbative 
QCD [4J. Our numerical results are presented in Figures I2ll3f5| and [6] for an illustrative 
purpose. As the results show, the real part of chiral magnetic conductivity stays to the 
value at u = for small u, contrary to the result in weak-coupling where it drops to 
| as soon as uj ^ [4J. As was already pointed out in ref. [3], it might be due to the 
strong interactions of the charge carriers. Consequently, the current response and charge 
asymmetry in RHIC plasma would be bigger than the weak-coupling result in ref. [3J. As 
a future direction, it might be worthwhile to generalize our calculations to the cases with 
non-zero electromagnetic charge density. 
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A Appendix : An easy method of performing numer- 
ical jobs 

We explain our way of solving the equation (I4.58P 

d r (-iurfo + rV(r)9 r / ) - iwr (d r f ) = , (A.105) 
d r (-iurg 1 + r 3 V(r)d r g 1 )-icjr(d r g 1 ) + t ^^f = , (A.106) 

to compute o~(u>) numerically. One first notices that it is straightforward to solve the 
equation ( 1A.105I) for fo numerically without any shooting ambiguity. For regularity at 
r = th where V(r) = 0, one simply considers r — > r# limit of (1A.105I) to deduce the 
regularity condition 

{drh){r H )={ 31777 T 9- W") > ( A ' 107 ) 

\r 6 H V'{r H ) - 2iur H ) 

and since a normalization of /o is not important, we can start at r = r# with /o( r /f) — 1 
and (d r fo) (ru) given in the above to solve (1A. 105ft numerically for r > r^. 

Once /o is found, there is a nice method of solving ( 1A.106I) for gi without further 
numerics. The idea is to transform the differential operator acting on g% into an integrable 
form, 

d r (-iur-+r 3 V(r)d r -)-iur(d r -)=P(r)d r (^R(r)d r (S(r)-)^ , (A.108) 

and comparing the both sides gives us the equations for the yet unknown functions P(r), 
R(r), and S(r) as 

PRS = r 3 V(r) , (A.109) 
p(R(d r S) + d r (RS) S j = d r (rV(r)) - 2iur , (A.llO) 

Pd r [R(d r S) J = -iu . (A.lll) 



Replacing R in the second equation (IA.110I) by using flA.109j) . one gets a nice simplification 



* ta i|) = -^o ' (A ' 112) 
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whose integration gives 

p-\r) = S-\r)e J - . (A.113) 

Using this and (1A.109I) to replace P and R in the last equation flA.lllj) . one gets a 
second order differential equation for S~ l , which turns out to be precisely the original 
homogeneous equation with the differential operator (1A.108j) . whose regular solution we 
already obtained as fo numerically. Therefore, one can simply let S^ 1 = fo and one finally 
has 

S-\r) = / (r) , 

n, f r dr' 

P -1 (r) = f (r)e J »(^WF) 

R(r) = r 3 V(r) (fo(r)) 2 . (A. 114) 

The equation (1A.106|) for g 1 subsequently takes a form 

P(r)d r (R(r)d r (S(r) gi )\=- l -^-fo , (A.115) 

whose first integration results in 

R(r)d r (S(r) 9l ) = -U2kQ j^dr* J^p^ , (A.116) 

where we fix an integration constant by considering r = r# where the left-hand side 
vanishes by R{th) = 0. The subsequent integration gives 




(A.117) 



where the second integration constant is fixed to make g\ normalizable at r — > oo. With 
the explicit choice of P, R and S in (1A.114I) . one obtains gi solely in terms of the known 
solution f as 

1i , r' dr" 

fr e ZW loo (r") a V(r") r r' (fn(r")) 2 ■ f" dr'" 

9l (r) = -iUnQMr) / dr' / dr" iM-i^e" 2 ^ ?^W> . 

Joo {r'YV{r'){fo{r')Y Jr H [r") 6 

(A.118) 

Having solved for gi, what one needs in order to find cr(u) is the near boundary 
expansion of <7i(r) in r — > oo. From V(oo) = 1 and the fact that /o(oo) is some finite 
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constant, it is indeed checked easily that g± has a behavior ~ -y as r — > oo, and more 
precisely 

Sl(r) ^f^r^(^ e --/r<^ 5 N )i + ... . (A.H9) 

V/o(oo) Jr H (r'y J r 2 

In conjunction with our formula for a (to) ( 14.591) and the result for the zero frequency limit 

one finally arrives at the useful expression 

which can be directly computed numerically once we find fo only. Our plots are generated 
using this technique. 

The above procedure can easily be repeated in the case of Sakai-Sigimoto model too, 
and we simply provide the results of each steps. For simplicity, we work in the unit 
2ty1 2 = 1. Once fo{U) is found numerically by requiring at the horizon U = Ut 

the differential operator acting on g\ can be transformed to an integrable form with 

S-\U) = f (U) , 

P-\U) = f {U)e - 2lu} ^ dU cm , 
R(U) = C(U)(f (U)) 2 e' 2iuJ ^ dU ' , (A.123) 
so that the equation for g\{U) is integrated as 

r U' JTT „ B{U") 

9i(U) = -fjfcfo(U) / dU' e - -J dU"F tu {U") (h{u")) 2 e^L «> ikA 

47r 2 C Joe C(U')(f (U')) Ju * V ' 

(A.124) 

where 

F « {U) -7ukw ■ (A ' 125) 

From the near boundary behavior C(U) — > ?/2 as C/ — > oo, one deduces without difficulty 
that 

- fl-#7-T f d^«^(/o(^))V^'- dU "^P>) , (A.126) 

\67r^C/o(oo) 7i/ T v 7 / [/a 
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from which in conjunction with fl 5 . 9 5 [) and (15 . 1 03j) . 



2 

a(cu = 0)=a = ^-N c , (A.127) 



one finally obtains the expression 



^1 = l — r d U> - a (f (U')) 2 e^S~ dU "^ , (A.128) 

solely in terms of fo(U)- 
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